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Even if a noninteracting system has zero Berry curvature everywhere in the Brillouin zone, it is possible to in¬ 
troduce interactions that stabilize a fractional Chem insulator. These interactions necessarily break time-reversal 
symmetry (either spontaneously or explicitly) and have the effect of altering the underlying band structure. We 
outline a number of ways in which this may be achieved, and show how similar interactions may also be used to 
create a (time-reversal symmetric) fractional topological insulator. While our approach is rigorous in the limit of 
long-range interactions, we show numerically that even for short-range interactions a fractional Chem insulator 
can be stabilized in a band with zero Berry curvature. 


I. INTRODUCTION 

There has been a great deal of recent interest in fractional 
Chern insulators (FCIs), which are the lattice analogs of frac¬ 
tional quantum Hall (FQH) states.® In addition to providing 
new experimental settings in which to explore topologically 
ordered phases of matter, FCIs demonstrate lattice effects that 
substantially change the underlying physics from the quantum 
Hall picture. 

The usual description of FCIs proceeds in the following 
way: Imagine starting with a so-called Chern band—a band 
which, when completely filled with noninteracting electrons 
would have nonzero quantized Hall conductivity. One next 
imagines partially filling the band and introducing an interac¬ 
tion which creates a gapped, many-body ground state. If the 
interaction is chosen correctly, the resulting state will have a 
fractionally quantized Hall conductivity. 

While this has proved to be a reasonable and successful 
approach to the problem, we believe it overemphasizes the 
structure of the noninteracting band, which may be irrelevant 
if the particles interact strongly. To make this point clear, we 
will construct an example that starts with a noninteracting sys¬ 
tem having zero Berry curvature everywhere in the Brillouin 
zone. Then, by adding an appropriately chosen interaction 
that breaks time reversal symmetry (either explicitly or spon¬ 
taneously), we will find that the ground state is an FCI. We 
will then show how similar interactions may be used to gener¬ 
ate a time reversal-symmetric, fractional quantum spin Hall 
insulator,® which is an example of a fractional topological 
insulator (FTI). We support our findings with numerical ev¬ 
idence from exact diagonalization studies. While our analytic 
arguments rely on the use of very long-range interactions, we 
find in our numerical work that even fairly short range interac¬ 
tions can form fractional Chern insulators in bands with zero 
Berry curvature. 


Although our construction is very general, let us restrict our 
attention to the case of spinless particles on a honeycomb. 
This is just a triangular lattice with two sites per unit cell, 
which we will call A and B sites. We consider a kinetic en¬ 
ergy, 

K = clc^,, ( 1 ) 

r,r' 

where r and r' are summed over the sites of the lattice and 
trr' describes the hopping amplitudes. In order to produce a 
Chern band, the hoppings (which are short-range) couple the 
sublattices, and phases are introduced which break time re¬ 
versal symmetry. Typically one tunes hoppings up to third 
nearest neighbors in order to get a very flat (although not per¬ 
fectly flat) Chern band. If desired, one can use further neigh¬ 
bor hoppings to make the bands extremely flat. For example, 
trr' may be chosen to describe the flattened Haldane honey¬ 
comb model.®® 

If one introduces a suitable short-range interaction (such as 
a nearest-neighbor repulsion), which we write as 

V = ^ ^ (2) 

rr' 

with fir = clcr being the particle number on site r, then the 
total Hamiltonian 

If = K + V (3) 

will produce an FCI ground state at the appropriate particle¬ 
filling fraction. For example, FCI states analogous to Laugh- 
lin states may be observed at zz = 1/3 for fermions and 
It = 1/2 for bosons. This conclusion has been established 
numerically by a number of groups.®^ 

III. FCI IN A BAND WITH ZERO BERRY CURVATURE 


II. ECI IN A BAND WITH NONZERO CHERN NUMBER 


We will first outline the conventional method for obtaining 
an FCI state from an underlying topological band structure. In 
Sec. Ill we will then attempt to reproduce this state, starting 
from a band with zero Berry curvature. 


Now consider a different band structure, described by a ki¬ 
netic energy K, which has zero Berry curvature throughout 
the Brillouin zone. We claim that we can switch on an inter¬ 
action, V, such that the ground state of this new model is also 
an FCI and is adiabatically connected to the FCI described 
in Sec. Our plan is to choose the interaction V such that 
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H = K + V mimics the Hamiltonian ([^, which is known to 
give an FCI. 

To be specific, let us write the noninteracting kinetic energy 
as 


We may now consider the interaction term 


17 = i ^ Trr' : Nrclc^, + clc^,Nr' 


(9) 


iT = ^irr' 4Cr') (4) 

r,r' 

where t^-r' describes new hopping parameters. For simplicity 
we will consider a hopping model which is diagonal in the 
sublattices, i.e., 

f if r, r' G A, 

U.r' = < ffr' if r, r' G B, (5) 

[ 0 otherwise. 

We also break the symmetry between sublattices by giving 
each one a different on-site energy In this way, we can 
arrange such that the two bands do not overlap in energy, the 
eigenstates are of the simple plane wave form and each 

band involves only one sublattice. Since the bands are com¬ 
pletely decoupled, for any translationally invariant hopping 
trr', the eigenfunctions have trivial structure and all bands 
have zero Berry curvature throughout the Brillouin zone (see 
Appendix [A| for the definition of this quantity). 

Now, certain interactions, which are short-range with a de¬ 
cay length ^ on the order of several lattice spacings, can be 
treated accurately within mean-field theory as long as ^ is 
much larger than the interparticle spacing. Consider an op¬ 
erator 

TVr = ^ V ./(r' - r)nr' (6) 

n ^' 

r' 

where n is the average density of electrons in the system. In 
this expression, / is some function (like a Gaussian) that is 
smooth and slowly decaying with length scale and has nor¬ 
malization 

^/(r'-r) = l. (7) 

r' 

By short-range with decay length we mean that the function 
f(r' — r) satisfies 

|/(r'-r)| < (8) 

for all sufficiently large |r' — r| and with finite constants C 
and 

In this way, the operator JVr measures the density of elec¬ 
trons in a region of length scale ^ and compares it to the av¬ 
erage density, n. If the state being operated on has approxi¬ 
mately uniform density, and if we take ^ to be large enough, it 
must be the case that Nr is extremely close to unity (since we 
are comparing the average density in a very large region to the 
real average density). We can then accurately approximate Nr 
by unity and do not need to treat it as an operator. As we take ^ 
larger and larger, this approximation becomes more and more 
accurate (assuming there is no phase separation). The fluctua¬ 
tions of Nr around this average may be treated perturbatively 
(see Appendix [B). 


where Trr' is an interaction strength (to be described below) 
and the entire expression is normal ordered. This operator 
is Hermitian and includes four-fermion interaction operators, 
although it is not generally a density-density interaction. The 
interaction is short-range in that it involves particles that are 
separated by a maximum distance on the order of In the 
mean-field limit. Nr becomes unity and we obtain 

U -^Y^^rr'clCr,, (10) 

rr' 

which is simply a hopping term. Thus, let us choose 

Trr' = irr' - trr', (H) 

where trr' were defined in Eqs. Q and ([T]), respec¬ 

tively, so that 

K + U^K ( 12 ) 

at the mean-field level. In this way, our scheme uses an inter¬ 
action to modify the “effective” band structure of the under¬ 
lying (topologically trivial) model. Note that since trr' breaks 
time reversal, so too does Trr' 

One might worry that this density-dependent interaction is 
unstable to phase separation, and indeed, it is fairly simple to 
write down an expression for Trr' for which this occurs. To 
make sure that this is always disfavored, we can add a diagonal 
term to the interaction strength, writing 

Trr' —t Trr' + MSrr', (13) 

where the second term may be absorbed into the amplitude 
Trr'- When added to Eq. this diagonal term generates a 
density-density interaction between a particle at r and its sur¬ 
rounding neighborhood, suppressing density fluctuations on 
the scale of By adjusting the strength of M, we can en¬ 
sure that complete phase separation is never favored, but that 
the small density fluctuations required for an ECI state are al¬ 
lowed. 

As mentioned previously, our mean-field picture is exact in 
the large ^ limit. Another concern is whether this mean-field 
state remains stable in the presence of fluctuations about this 
limit. In Appendix [B] we discuss the perturbative corrections 
due to the interaction term U and justify the stability of the 
ECI state. We then provide numerical evidence in support of 
the ECI state in Sec.lVl 

The interaction U changes the effective band structure of 
the model, so if we now use the total interaction 

V = V + U, (14) 

then the Hamiltonian H = K + V completely mimics the 
above H [Eq. ([^] and will have an ECI ground state. As re¬ 
quired, the noninteracting band structure of K has zero Berry 
curvature and the topological properties are induced purely by 
the interaction term, V. 
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The idea underpinning this approach is that hopping terms 
are essentially a subset of all possible interaction terms. By in¬ 
troducing an auxiliary operator, in this case a local density op¬ 
erator, that acts as a number on the ground state, an interaction 
can be written down that modifies the effective single-particle 
band structure. The trade-off for this is that the interactions 
are no longer density-density, and have a range that is finite 
but may cover several unit-cell lengths. 

In the expression above we have broken time-reversal sym¬ 
metry by hand with our choice of interaction U. If preferred, it 
is possible to generate an FCI state spontaneously with a time- 
reversal-invariant interaction, as discussed in the next section. 
This requires us to introduce a spin-1/2 degree of freedom but 
will also allow quantum spin Hall insulator (or FTI) states to 
be realized. 


IV. TIME-REVERSAL-INVARIANT INTERACTIONS 
A. Spin-dependent band models 

One objection to the approach outlined above is that it 
breaks time-reversal symmetry explicitly through the interac¬ 
tion U. In order to use a time-reversal-symmetric interaction, 
we now promote the spinless fermions to spin-1/2 fermions. 
These obey the usual interrelations under time reversal, which 
can be written as 


0Cra0 ^=iTl^,Cra', (15) 

where is the Pauli y-matrix, a and a' are spin labels, and 
0 is the time reversal operator.!^ We have also included a spin 
index on the field operators Cra- 

The underlying topologically trivial bands must now have a 
spin dependence. For simplicity, we will treat spin as a good 
quantum number and write the kinetic-energy term as 

= (16) 


The framework outlined above leads to a Hamiltonian that 
conserves spin. For example, K' could describe the kinetic 
ener^ of the Kane-Mele model without the Rashba spin-orbit 
term.l^lln principle, our approach could be generalized to time- 
reversal invariant models that do not conserve spin, such as the 
full Kane-Mele model. However, in that case, the interactions 
we describe below would involve complicated terms that mix 
both physical spin species; we will restrict our discussion to 
the simpler, spin-conserving case. 

To generate a fractional state, we would partially fill the 
lowest band for each spin species and turn on a short-range 
interaction, 

V — Uj.j./^7Zr(T^r'(T- (19) 

rr'(7 

This is diagonal in the spin index, so at this stage there is no 
interaction between spins. 

The total Hamiltonian 

H' = K' + V' (20) 


would generate a fractional topological state, but the overall 
Chern and spin Chern numbers would depend on the filling 
fraction of each of the two lowest spin bands. If both bands 
are filled equally, say at i/ = 1/3 each, then the system will 
retain its time-reversal symmetry and will form a fractional 
topological insulator state. If one band remains empty but 
the other is (fractionally) filled, then the system could sponta¬ 
neously break time-reversal symmetry and form an FCI state, 
depending on the interaction V. The relative stability of these 


two possibilities will be discussed in Sec. IV C 


B. Topological phases in bands with zero Berry curvature 


We now follow a similar approach to Sec. Ill to write down 
an interaction which changes the effective band structure from 
K' to K'. We use the local-density operator Nj. as defined in 
Eq. (§, but the lattice site density operator now sums over 
both spin species. 


although in general we could have included spin-flipping 
terms. As before, we will assume that the different sublattices 
have different on-site energies. Then, following the arguments 


of Sec. Ill each of the four bands will have exactly zero Berry 
curvature. 

For the hopping parameters of our target Chern band model, 
we take /rr' from Eq. Q for one spin species but, crucially, 
invert the Chern number for the other spin species by taking 
the conjugate hopping. In other words, we take 


f' — f / 

~ prr'] 


(17) 

(18) 


and set any spin-flipping terms to zero. This defines a new 
topologically nontrivial kinetic energy, K', which one can 
easily show is symmetric under time reversal symmetry. This 
kinetic energy therefore generates a low-energy band for each 
spin species, and these have equal and opposite (nonzero) 
Chern numbers when filled. 


fir — ^ ^ ffrc- (21) 

<7 

With this, we write 

u'=Y. 

rr' ( 7 ( 7 ' 

where any diagonal interaction-hopping terms necessary to 
prevent phase separation have been absorbed into T^^,. 
Choosing 

T'r' = i'rr' " Kr' ( 23 ) 


N, 


r^ra^r' a‘ 


+ cl 


.CrV'^r' 


( 22 ) 


and following our previous reasoning, in the mean-field limit 
one finds 


K' + U' ^ K', 


(24) 


which generates a spin-Chern band for each spin species. Fi¬ 
nally, using the complete interaction U' -f V', it is possible to 
generate a fractional topological state from the noninteracting, 
zero-Berry-curvature system. 
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C. Tuning ferromagnetism 


tened Chem band hopping K is defined by 


We noted previously that Hamiltonian ( |20l l could gener¬ 
ate either an FCI or an FTI state, depending on whether or 
not time reversal symmetry was spontaneously broken. For 
example, if the lowest two bands were filled with density 
(i/^, v^) = (1/3,1/3), then a state with zero Chern number 
but nonzero spin Chem number would be produced. On the 
other hand, if the fillings were = (1/3,0), then an 

FCI state with nonzero Chern number would be produced.^! 

These states occur at different overall fillings but also de¬ 
pend on whether the system is susceptible to spontaneous fer¬ 
romagnetism. In order to tune this susceptibility, we introduce 
the spin-spin interaction 

l"spin = a^(^rT-iVr;)', (25) 

r 

where the operators 

Nya = 3 V] /(r - r')hya (26) 

r' 

calculate the local spin density in a region of linear size 
For a > 0, differences in spin density will lead to an energy 
penalty, encouraging the fermions to be split equally between 
both spin species. For a < 0, spontaneous ferromagnetism 
will be energetically favored, and fermion spins will tend to 
align. By tuning a, it is therefore possible to produce both FTI 
and FCI states from a topologically trivial system by adding a 
time-reversal-invariant interaction. 

The philosophy of this section is essentially the same as the 
spinless fermion case; by choosing a suitable interaction, the 
effective band structure of the single-particle problem can be 
altered. Using spin-1/2 fermions allows us to do this without 
breaking symmetry explicitly, and we can thus generate both 
FCI and FTI phases. 


V. NUMERICAL RESULTS 


To support the arguments above, we now give some finite- 
size numerical results which show an FCI state that has been 
generated using interactions of the form described above. Al¬ 
though our above arguments are rigorous in the limit of very 
long-range interactions, we will find that even for a short- 
range interaction, we can stabilize a fractional Chern insulator 
in a band with zero Berry curvature. 

We will consider a model of bosons hopping on a hon¬ 
eycomb, forming bands with zero Berry curvature (as in 
Sec. mi- We will consider a case where the interaction 
strength (which determines the scale of K) is much greater 
than the hopping bandwidth K. For simplicity we can set 
tyy' — K = 0, which is appropriate in this limit. For the 
kinetic energy K that we wish to simulate, we consider the 
flattened Haldane model, which is outlined in Ref.|5] The flat- 


Khh = —ii ^ \c\,c^ -\- H.c. 

(rr') 

-h X! [ci'Cr + H.c. 

(((rr'»> 


t2 ^ -f H.C. 

((rr'» 

+ M^4cr, (27) 


where ti = l,f 2 = 0.6, fa = —0.58, (j) = OAtt, and HH 
stands for Haldane honeycomb. This has a flatness ratio of 
about 50 and exhibits av = 1/2 FCI state that remains stable 
even in the presence of first- and second-neighbor density- 
density interactions.!^ In order to prevent phase separation, we 
have also added a diagonal term with strength M, which in 
the noninteracting model appears as a chemical potential. 

For simplicity, we assume we have hardcore bosons, and 
we take the function /(r) to be a simple top hat. 


/(r' - r) 


n(l) 

0 


if |r' — r| < Z, 
otherwise. 


(28) 


In this expression, I is a length scale that covers a small num¬ 
ber of nearest neighbors, and n{l) is the number of sites that 
lie within the top hat, so that the function is normalized. Us¬ 
ing Eq. (|^, we can now write the flattened Haldane model as 
the interaction 


Ukk = 2 


HH 


Nyclc^i + clc^iNy 


(29) 


where T™ contains all the hopping parameters from Eq. ( [ZT] ) 
(including the diagonal M term). Since we are considering 
hardcore bosons, we do not need to add any additional inter¬ 
actions in order to stabilize the ECI state at half filling. 

We carry out exact diagonalization of Uhh on a 24-site lat¬ 
tice with 6 particles and on a 32-site lattice with 8 particles. 
In both cases we choose the size of the top hat, defined by 
I, to extend up to third-nearest neighbors; note that this is no 
farther than the hopping distance. Eor the 24-site case we set 
M — 2, and for the 32-site case we set M = 3, chosen to give 
the largest many-body gap to ground-state splitting ratio. In 
both cases, we find strong numerical evidence to suggest that 
the ground state is a ;/ = 1/2 ECI. 

In Eig. [T]we plot the low-lying energy levels for each lat¬ 
tice size as magnetic flux (j) is threaded through one handle of 
the torus. In each case we find a twofold-degenerate ground 
state that is separated from excited states by a large many- 
body gap and that occurs at the momentum values predicted 
for a:/ = 1/2 ECI by the generalized Pauli principle.!^ As flux 
is threaded through the handles of the torus, the two ground 
states for the 24-site case evolve into each other with a level 
crossing. Eor the 32-site case, the two ground states are in 
the same momentum sector, and display an avoided crossing 
under a flux insertion. 

We have also calculated the particle entan gleme nt spec¬ 
trum for various partitions for each lattice sizeP^El In Eig. 
we show the particle entanglement spectrum for the 32-site 
lattice, tracing out three of the eight particles (we form the 
density matrix from an incoherent superposition of the two 
degenerate ground states, p = 5 X) I'*/*) motivated 
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0/(2i) i>/(2;r) 

Figure 1. (Color online). Low-energy spectra as flux (f) is inserted 
through one handle of the torus for the half-filled honeycomb lattice 
with (a) 24 sites and (b) 32 sites. Insets show a close-up of the ground 
state energy splitting. This figure gives evidence that the ground state 
of our model Hamiltonian [bosons on a honeycomb with zero hop¬ 
ping and interaction given by Eq. |29| is a fractional Chem insulator. 
See main text for details of the model. 


Ni = N 2 = 4, Np = 8, = 5, Np = 3 



N2ki+k2 

Figure 2. (Color online). Particle entanglement spectrum for N = 8 
hardcore bosons on a 32-site honeycomb, with Nb = 3 particles 
traced out. Ni and N 2 are the numbers of unit cells in each direc¬ 
tion, with fei and k 2 being the corresponding momentum integers 
in the Brillouin zone. The number of states below the dashed line 
matches the expected quasihole counting rules (22 states in each sec¬ 
tor, giving a total of 352 states.). This figure gives evidence that the 
ground state of our model Hamiltonian [bosons on a honeycomb with 
zero hopping and interaction given by Eq. is a fractional Chem 
insulator. See main text for details of the model. 


in Ref. [Ull. The entanglement spectrum shows a clear en¬ 
tanglement gap, and the number of states below this gap 
obeys the expected (1,2)-admissible quasihole counting rules 
for each momentum sector.E^l This strongly suggests that our 
interaction-stabilized state is in the same universality class as 
other v = 1/2 FQH and FCI states. 


VI. DISCUSSION 

In this paper we have proposed a type of interaction that 
may be used to generate fractional topological phases from 
a topologically trivial band structure. The underlying idea 
throughout is that hopping, i.e. kinetic-energy terms, are ef¬ 
fectively a subset of the possible interaction terms. As such, it 
is possible to change the single-particle spectrum using an in¬ 
teraction rather than a kinetic-energy term. We have outlined 
how this may be achieved by interacting with a local density 
in a region of size it is possible that other approaches exist 
which use the same philosophy. 

Although we have used a formally short-range expression 
for Nr in the arguments above, our approach should also ap¬ 
ply more generally (including, for example, cases where 
decays algebraically). The essential requirement is that the 
function /(r' — r) be normalized as in Eq. Q, although it will 
also need to be suitably well behaved so that it does not favor 
charge-density-wave states. El 

We note, however, that to change the Chern number of 
a band using an interaction, the interaction must be strong 
enough to mix the single-particle bands; that is, it must be on 
the order of the gap between the bands before the interaction 
is added. It is often the case that FCIs are studied for flat bands 
in the very weak interaction limit. For this case our strategy 
will not work since a very weak interaction would not be able 
to change the effective Chern number of the band. However, 
one could also consider a more complicated strategy where the 
lattice translational symmetry is broken down further, enlarg¬ 
ing the unit cell such that a band with zero Chern number is 
broken into two bands with opposing nonzero Chern numbers. 
This could be done with an arbitrarily weak perturbation. 

Although the interactions we consider are not of the 
density-density type, density-dependent hopping interactions 
are predicted to arise in systems of cold atoms under cer¬ 
tain conditions. They may arise as higher-order terms in cer¬ 
tain tight-binding approximation^ and are also suggested 
to occur in e ffective Hamiltonians when band mixing is 
important.EnSI More directly, density-assisted hopping terms 
may be induced using a periodic ally d riven magnetic field 
to modulate interaction strengths, EEH an approach which 
may also be used to produce density-dependent synthetic 
gauge fields.El In addition, vacancy-assisted hoppings have 
recently been discussed in the context of kinetically con¬ 
strained models, where they may also lead to topologically 
ordered phases.®! 

Nonetheless, these proposals have yet to be realized in a 
laboratory, and in their current formulation would not be able 
to realize the specific local density Hamiltonians discussed in 
this paper. For the time being, experimental realizations of 
FCI states seem most likely to be achieved through the con¬ 
ventional Chern band route, a direction in which substantial 
progress has already been made.®3®l 

Finally, it is interesting to note that our construction pro¬ 
vides a nice counterexample to the claims of Ref. (which 
have also been disproven in another manner by the current au¬ 
thors in Ref. 1251) . Reference |3] claims to calculate the Hall 
conductivity of an FCI as an integral of the Berry curvature of 
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the noninteracting bands multiplied by some occupancies over 
the Brillouin zone. It is clear that this cannot be correct since 
in our construction the Berry curvature of the noninteracting 
bands is strictly zero (so that Ref. |31 would always predict 
zero Hall conductivity) whereas, due to the (time-reversal- 
breaking) interaction, the system is an FCI with quantized, 
nonzero. Hall conductivity. 


of the underlying torus. If the many-body ground state is writ¬ 
ten as Irk) and we thread flux (pi and (/)2 through each handle, 
then the many-body Berry curvature is defined through 


F{(j)l,(t)2) 




\ 

d(j)2 dpi) 


■ (A4) 


The many-body Chern number can then be calculated by inte¬ 
grating the above quantity over pi and p 2 . 
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Appendix A: Berry Curvature and Berry Matrix 


Appendix B: Stability Beyond Mean-Eield Theory 

In Sec.|I^we introduced an interaction that reproduces the 
kinetic energy of a target Chern band model exactly in the 
mean-field limit. In this appendix we will consider perturba¬ 
tive corrections to this mean-field picture and justify the sta¬ 
bility of the PCI ground state. 

When the correlation length ^ tends to infinity (or reaches 
the system size in the finite-size case), the “local” density op¬ 
erator Nr can be replaced with unity. In general, the operator 
U can be written as U = Umf + SU, with 


In this appendix we will briefly outline the definitions of 
Berry curvature that we use in the main text. Pirst, we write 
the single-particle Bloch solution for the ath band as 

pJakir) = e**"'’uak(r), (Al) 


where Mak(r + R) = Uak(r) is periodic under lattice transla¬ 
tions. Por a single band, the (gauge-invariant) Berry curvature 
is defined through 


F“(k) 



dUak \ 
/’ 


(A2) 


where, for a two-dimensional system in the xy plane, we mean 
the kz component of the cross product. Integrating this quan¬ 
tity over the Brillouin zone (and dividing by 27r) gives the 
Chern number, which is proportional to the transverse Hall 
conductivity of the filled band and must take an integer value. 

When there are several degenerate (or nearly degenerate) 
bands, the generalization of the Berry curvature is the Berry 
matrix. 



where the sum over c includes all bands in the degenerate mul- 
tiplet. The Berry matrix is not gauge invariant itself, but mea¬ 
surable quantities are contained within gauge-invariant ex¬ 
pressions like the trace. We note that the components of this 
matrix vanish trivially if the multiplet contains all bands of the 
system. 

In an PCI state, the noninteracting band picture breaks 
down and we can no longer calculate the Chern number from 
the single-particle properties. Instead, the many-body Berry 
curvature is calculated by threading flux through the handles 


t^MF = ^ Tr, 


'C'C,, 


SU = : SNrclc^, +clCj.,SNr' 


(Bl) 

(B2) 


and 

SNr = iVr - 1 = ^ V f{r' - r) [hr - h]. (B3) 

n 

r’ 


In this way, the perturbation SU depends on the density fluctu¬ 
ations in a region of linear dimension scaled by the average 
particle density h. 

In order to estimate the size of the corrections due to <5(7, 
we first recall that for a general, uncorrelated Poisson process, 
density fluctuations in an area of linear size ^ scale as 8Nr ^ 
1 /^. Quantum Hall states are more correlated than this, so we 
expect their density to fluctuations to scale as 8Nr 1/C 
most. Assuming that this property transfers to systems on a 
lattice, we expect SNr acting on an PCI ground state to have 
an effect that is no larger than 0(1/C). 

The action of the hopping component of SU is bounded 
by the maximum bandwidth W, multiplied by the number 
of neighboring sites and orbitals involved in the sum h. At 
pth order in perturbation theory, the action of SU on the PCI 
ground state is then bounded by 


su^p^ < 



P 


(B4) 


with A being a constant that depends on energy denominators 
and scale factors of order unity. By choosing C large enough, 
we can ensure that perturbative corrections of this form are 
negligible and convergent. 

However, our assumption that SNr ~ 1/C is valid only at 
low orders in perturbation theory when the state being acted 
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upon is close to a (homogeneous) FCI ground state. After ap¬ 
plying SU many times, perturbative hopping processes may 
generate virtual phase-separated states, which have 6Nr ~ 1. 
Assuming that the hopping described by T^r' can transfer par¬ 
ticles by up to L sites, virtual states that move particles by 
linear dimension ^ will first start to appear at order ^/Lin per¬ 
turbation theory. A finite density of particles over a region of 
area may then be moved, hence significantly changing the 
value of N at order in perturbation theory, where 


C is some constant. We can push these contributions to arbi¬ 
trarily high order in perturbation theory by making ^ larger. 
Additionally, we can suppress these perturbative corrections 
with a large energy denominator by choosing M, the diagonal 
interaction term in Eq. ([T3J, to be large and by adding addi¬ 
tional density-density interactions between nearby neighbors. 
Finally, we note that our numerical work confirms that even 
for small fractional Chem insulators may be stabilized. 
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